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The complex interactions between internal motor pressure oscillations resulting from 
vortex shedding, the motor’s internal acoustic modes, and the motor’s structural vibration 
modes were assessed for the Space Shuttle four-segment booster Reusable Solid Rocket 
Motor and for the five-segment engineering test motor ETM-3. Two approaches were 
applied: 1) a predictive procedure based on numerically solving modal representations of a 
solid rocket motor’s acoustic equations of motion and 2) a computational fluid dynamics two- 
dimensional axi-symmetric large eddy simulation at discrete motor burn times. 

Nomenclature 


A = aft segment 
AC = aft center segment 
f = frequency (Hz) 

F = forward segment 
FC = forward center segment 

i = segment index; also used as a location or spatial index 

j = vortex stage number index 

L = length 

m = mode index 

p = pressure 

S = Strouhal number 

t = time and its derivatives d/dt and dVdt^ 

U = axial gas velocity 

to - frequency (radians/second) 

C, = viscous damping 
cp = mode shape coefficient 
x = distance 


I. Introduction 

Large segmented solid rocket motors have a history of exhibiting pressure oscillations with corresponding thrust 
oscillations. The Space Shuttle reusable solid rocket motor (RSRM) also exhibits low level pressure and thrust 
oscillations. Pressure and thrust oscillations were noted during the first ground test 1 and have been monitored and 
studied since then. 

Typical data from a recent RSRM static test is shown in Figs. 1 and 2. Measured head-end pressure typically 
contains quasi-sinusoidal signals associated with the internal longitudinal acoustic or organ pipe resonances of the 
motor (denoted as 1-L, 2-L, etc.). The 1-L and 2-L frequencies are near 15 and 29 Hz, respectively, and vary only 
slightly during motor bum. These frequencies are shown in Fig. 1 . Typically these oscillations start at about 40 
seconds into motor bum. 
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Figure 2: Measured Axial Thrust Data from RSRM FSM-8 Static Test 


Pressure oscillations in large, segmented, solid rocket motors are often the result of a vortex shedding 
phenomenon. Shearing action of the exhaust gases as the flow passes across internal segment slots and protruding 
inhibitor stubs causes vortices to form. Flatau presented the hypothesis that the shearing flow gives rise to vortices, 
which travel downstream and interact with a downstream baffle such as a protruding inhibitor stub. 2 ' 1 This action 
creates a sound or whistling noise. The frequency of the vortex shedding sound is a function of the axial gas velocity 
and the distance between inhibitors. 

fjj = Sjj * Uj / Lj (1) 

where 

i = segment number 

j = ‘"stage number,” which some have suggested is the number of vortices that fit in the space between 
upstream and downstream baffles, j has some characteristics of a harmonic. 

Vortex shedding frequencies decrease with motor bum time, since the axial gas velocity decreases with time as 
the propellant grains inner diameter increases. 

Vortex shedding excitations have the potential to interact with the acoustic and structural modes of the motor 
(Fig. 3). 
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Two approaches were 
applied to die problem of 
predicting and understanding 
pressure and thrust oscillations 
in a solid rocket motor. The first 
approach was a time domain 
simulation that numerically 
solved modal representations of 
a solid rocket motor’s acoustic 
equations of motion. This 
simulation was anchored using 
an extensive database of four- 
segment RSRM measured data. 

The simulation was then 
extended to the five-segment 
engineering test motor (ETM-3). 

The second approach was a Figure 3: RSRM Acoustic, Structural, and Vortex Shedding Interactions 
computational fluid dynamics 

(CFD) two-dimensional (2-D) axi-symmetric large eddy simulation (LES) using the turbulence model option in 
FLUENT®. CFD models were developed at discrete motor bum times. Initially, the CFD effort was intended to 
better define some of the weakly defined parameters of the first approach described above. However, the CFD 
modeling effort showed great promise in fundamentally understanding and predicting the internal flow behavior of 
the RSRM, ETM-3, and other solid rocket motors. 
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IL Simulating Vortex Shedding in the Time Domain 

The first approach was the development of a predictive procedure programmed on PV-Wave® software and 
based on numerically solving modal representations of a solid rocket motor’s acoustic equations of motion at time 
intervals for the entire duration of motor bum. The first three longitudinal acoustic modes of the motors were 
implemented in the solution. Acoustic mode data for selected motor bum times for both motors were obtained from 
both the Standard Stability Prediction Program (SPP) and NASTRAN acoustic models. 

Excitation to the system was assumed to be from sinusoidal sound (vortex shedding) sources located at each 
segment propellant slot and from distributed random (combustion noise) excitations along the length of the motor. 
The sinusoidal sound sources and their harmonics vary in both frequency and amplitude with motor bum time, and 
were derived from motor ballistic predictions. Strouhal number values for the harmonics were based on sub-scale 
cold flow simulations of the RSRM conducted by Flatau. 2 The analytical procedure and model parameters such as 
sound source strengths, relative harmonic strengths, and acoustic damping were anchored using RSRM static test 
data. 

The resulting simulated pressure oscillation time histories were then applied to NASTRAN structural models of 
the motor in the static test stand for discrete motor bum times. The resulting motor dynamic thrust and acceleration 
responses were calculated using Fast Fourier Transform (FFT) techniques to generate waterfall and contour plots of 
the simulated pressure, thrust, and acceleration time histories. The structural models were validated with pre-test and 
post-test modal surveys of the motor in the static test facility. Good correlation to measured static test data was 
obtained. 

This new predictive technique was used to make worst-case, pre-test dynamic pressure and thrust oscillation 
predictions for the five-segment ETM-3 static test. This process bounded pre-test loads of the motor and test facility 
compared to structural capabilities. 

The basic analytical approach was to solve the one-dimensional (1-D) partial differential equation of the acoustic 
system by assuming a separation of variables solution. 


Pi(x,t) = p i (x)p j (t) 


The dynamic response is the sum of the modal responses. 


P 1 ( x ' t ) = S m (p,) m (x) (Pi) m (t) 


( 2 ) 

( 3 ) 


4 

American Institute of Aeronautics and Astronautics 



Dynamic pressure excitation at the ith location was assumed to be of the form 


Pi(t) = Py sin(27t fjj t) (4) 

where 

f {j = Sy Uj / Lj and j denotes the stage number or the j A harmonic. (5) 

The head-end (x^HE) modal pressure (p HE ) m (t) resulting from the ith location excitation is obtained by numerically 
integrating the second-order single degree-of- freedom modal equation 

M& 2 (PmvXn + 2 ^ (0 m d/dt (pHE iXn + 0) m 2 (Png i)ffi = Pi (t) (6) 

This yields 

The mode shape ratio between the itb location and the head-end is (p^ i ) m 

(Pi)m( X ) = (pHE,i)m = (^HE ^ ^iXn (2) 

Therefore, 

(«) 

For multiple excitation locations, linear superposition was used 

PH E (t) = S i P HE ,(t) (9) 

Random excitations were applied to simulate the combustion noise effects seen in test data 

= ( 10 ) 


This procedure was programmed using PV-Wave® software. The resulting p HE time history was processed 
as if it were measured data. 

Numerous assumptions and weakly defined parameters are involved in this simulation. The RSRM pressure 
oscillation observed behavior was used to refine the assumptions and define the various simulation parameters. 

In the vortex shedding frequency equation Eq. (1), Sj was assumed to = 0.85 * j for j= 1-4 (j>4 were not observed 
in RSRM test data). Flatau also used j = 1-4 with S T = 0.9, S 2 = 1.6, S 3 = 2.5, S 4 = 3.5. 2 U { are 1-D port velocities at 
the slots obtained from 1-D ballistic runs (Fig. 4). 
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Figure 4: RSRM Velocity vs. Time 
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are distances between the slots (or aft slot to nozzle length for the aft segment). For RSRM: = L 2 = L3 = 320 

in., L 4 = 330 in. 

Vortex sound locations were assumed to be at the downstream slots (or nozzle for aft segment) consistent with 
the hypothesis of hole-tone generated sound. The RSRM has three such locations: 1 = FC-AC, 2 = AC- A, and 
3 = A-N (Fig. 5). 



Figure 5: RSRM Vortex Sound Locations 

The vortex sound pressure strength at each of the i* slots is assumed to be (p s )i = a, (1/2 p, Ui 2 ), where ai was 
obtained by matching RSRM data, a 2 ~ ai (U 2 )/Ui) , and a 3 = aj (U 3 )/Ui). This eliminates arbitrarily selecting three 
parameters. From RSRM data, it was noticed that the apparent vortex sound pressure strengths were approximately 
related by the above velocity ratios. Weighting factors on stage numbers and on time were applied as suggested by 
RSRM measured data. 

The random sound pressure locations were distributed along the length of the motor. For the i* excitation location, 
(pr)i = [ Pi (1/2 Pi Ui 2 ) + p 2 Aj ], where pi and p 2 were obtained by matching RSRM data and (1/2 p; Ui 2 ) = mid-segment 
dynamic pressure and A* = mid-segment surface area. Ui and Ai vary with time. The two factors comprising pt are an 
attempt to simulate the random sound originating both from aerodynamic flow (decreases with decreasing port velocity) and 
from propellant combustion (increases/decreases with increasing/decreasing bum area). A random time-weighting 
function was also applied as suggested by RSRM measured data. 

The longitudinal internal acoustic modes of the RSRM have been extensively studied throughout the program. 
Figure 6 summarizes some of this history. 
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Figure 6: RSRM Longitudinal Acoustic Modal Characteristics 


For the RSRM simulation, acoustic natural frequencies versus bum time were obtained from a three-dimensional 
(3-D) acoustic model followed by curve-fitting versus bum time to generate a continuous function. 

RSRM mode shapes versus time were obtained from 1-D SPP runs. Mode shape coefficients at discrete bum 
times were curve-fit to create a continuous mode shape variation with bum time. 

Acoustic damping was also assumed as a function of time and reflects a variation with gas density as motor 
pressure changes. Damping values were anchored to RSRM data. 

The motor in its test stand has numerous natural frequencies in the frequency range of the acoustic modes. Most 
are motor bending and/or torsional in nature and also involve circumferential wave or case shell motion. Several 
structural modes involve significant axial motion and are susceptible to being excited by longitudinal acoustic 
excitations. Numerous analytical models of the motor in its test stand have been developed over the years. Some of 
these models are finite element models of very high fidelity and represent the motor at several discrete bum times 
(Fig. 7). In addition, several modal surveys have been conducted to provide validation of analytical structural 
models and to provide modal damping information. 



Figure 7: Structural Finite Element Model of RSRM in T-97 Test Stand 
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The governing equations of motion were programmed and numerically integrated using PV-Wave® software. A 
time history of head-end pressure was generated at 0.001 -second time intervals. The resulting data were processed 
as if it were actual pressure data. The results are shown in Fig. 8; compare to Fig. 1. Although the amplitude scales 
are slightly different, the simulation is similar to the measured data. 
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b) Waterfall Plot: Simulated Dynamic Signal 



c) Contour Plot: Simulated Dynamic Signal 


Figure 8: Simulated RSRM Head-End Pressure 
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This same simulation procedure was applied to the five-segment ETM-3 motor as a predictive tool to bound 
pressure and thrust oscillation amplitudes, frequencies, and occurrence times. Vortex shedding frequencies were 
adversely shifted to create several worst-case conditions. Figures 10 and 11 present worst-case 1-L pressure and 
thrust simulations. 



ETU-3 Simulation (15July03> 



Frequency (Hz) 


a) Time History: Simulated Dynamic Signal b) Waterfall Plot: Simulated Dynamic Signal 



0 10 20 30 40 SC 


Frequency (Hz) 

c) Contour Plot: Simulated Dynamic Signal 


Figure 10: Simulated ETM-3 Worst-Case 1-L Head-End Pressure 
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Test Stand Thrust (Mlbf) 


ETM-3 MEOP Utah Thrust Plus Worst-Case 1L Dynamic Interaction 
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b) Waterfall Plot: Simulated Dynamic Signal 


c) Contour Plot: Simulated Dynamic Signal 


Figure 11: Simulated ETM-3 Axial Thrust for Worst-Case 1-L Pressure 
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III. Computational Fluid Dynamics Simulations 


The second approach was to use the 2-D axi-symmetric LES turbulence model option in FLUENT® at three 
different motor bum times. Although, the model does not allow for 3-D vortex stretching, it does allow the entire 
motor to be efficiently modeled with all the inhibitor stubs, with a reasonable number of grid cells and in a 
reasonable amount of computer time. The inhibitor stubs were modeled using the best information available for the 
height, thickness, and bending geometry. 

Noting that a 3-D LES model had been successfully used for the French Cold Flow Test, a 2-D planar LES 
model using FLUENT® was constructed for the same cold flow simulation. 4, 5 ’ 6 Results of the FLUENT® model 
were comparable to the previous work. 

The LES turbulence model is an unsteady flow model using a small time-step on the order of 5E-05 seconds. The 
size of the model is on the order of the size of a steady-state Reynolds Averaged Navier-Stokes (RANS) model or 
about 80,000 cells. Within the 2-D axi-symmetric option, the LES model uses the turbulent standard wall function. 
The FLUENT® solution controls were set as second-order for pressure and density. The central differencing scheme 
was recommended for both the momentum and energy solver. However, it was found that the LES model results 
looked reasonable with the power-law for the energy solver and it ran faster with better convergence. The sub-grid- 
scale model used the Smagorinsky- Lilly option. 

The steady- state RANS model included the propellant bum-rate power-law as a function of local pressure. 
However, this model was turned off when the LES model was turned on. Otherwise, the LES model took too long to 
converge. Even though, the bum-rate power-law is absent in the LES model, some of the effects are partially 
accounted for, because the mass addition does vary down the bore approximately with the mean pressure drop. It 
was found that after the LES model was activated, the model would have a slow transient for about 1.5 seconds 
before settling down to a quasi-steady pattern in which pressure oscillations were occurring about a mean pressure. 
After the slow transient decay, the pressure monitor data could be used for the FFT and contour variables could be 
used for plotting internal motor descriptions. All of the FLUENT® runs used the single-phase model or a 
homogeneous gas without discrete particles. 

Figure 12 is a full view contour plot of the vorticity for the RSRM and ETM-3. Careful inspection shows that the 
vortices from the forward inhibitor, which is the longest, remain the most distinct as they move down the bore. 
Because of the length of the forward inhibitor, these vortices are produced close to the motor centerline. Since these 
vortices travel the full length of the bore, they have time to migrate closer to the centerline and tend to dissipate by 
the time they reach the nozzle nose. 



a) RSRM at 100 Seconds 



b) ETM-3 at 100 Seconds 
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Figure 12: RSRM and ETM-3 Vorticity Contour Plots 
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A close-up view of the RSRM and ETM-3 forward inhibitor stub is shown in Figs. 13 and 14, respectively. 
Inspection shows that the vortices from the forward inhibitor do not interact significantly with the wall instability 
layer. Movies made from these still plots clearly show that the vortices rapidly leaving the inhibitor stub tip tend to 
pair and revolve around each other until they coalesce into one larger but weaker vortex traveling downstream 
toward the nozzle. The revolution of two vortices with like rotation about a common origin is in keeping with the 
classical inviscid analytical solution of two vortex filaments with like spin. 7,8,9 



Figure 13: RSRM Forward Inhibitor Vorticity View at 100 Seconds 



Figure 14: ETM-3 Forward Inhibitor Vorticity View at 100 Seconds 


The general explanation for the phenomena of the sustained vortex shedding pressure oscillations is that large 
vortex structures interact with the nozzle nose causing a pressure wave to move upstream, which excites the 1-L 
acoustic mode of the motor. Flandro presented an early description of this process for solid rocket motors. 10 The 
general topic of flow structure impingement on solid surfaces is very wide ranging. Flandro also presented a good 
survey of this subject. 11 Another source of sound was treated analytically by Rockwell, which shows that vortices 
moving around obstacles can produce acoustic waves. 11 

Since the interaction of vortex structures with the nozzle is important, a sequence of close-up views of the RSRM 
center and aft inhibitor stubs are shown in Fig. 58. 

Each inset in Fig. 15 has been given the label of Frame No. 1 to Frame No. 6. Each frame is 0.012 seconds apart, 
so that six frames constitute about one cycle of the first harmonic. These are just some of the frames to make one 
movie. Each frame shown here actually represents every fourth frame in the movie. The movie goes for several 
cycles so that 1 09 frames are used. 

Figure 15 shows the very active center inhibitor stub is rapidly shedding vortices that interact with the 
downstream wall instability layer. Basically, there are three different types of vorticity regions listed as follows: 
1) the wall instability layer, 2) the vortices shedding from inhibitor stubs, and 3) vortices shed from flow over 
propellant segment cavities. All three types have also been documented in literature. 13, 14, 15 
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It is interesting to track the large vorticity region (magenta color) just upstream of the aft inhibitor and another 
large vorticity region midway between the aft inhibitor and the nozzle nose. The second region is just down stream 
of a small low vorticity (light blue) region. Proceeding to Frame No. 6 it can be seen that the first vorticity structure 
is stretched and pulled over the aft inhibitor stub, while the second structure is stretched and pulled towards the 
nozzle nose. The latter structure appears to either slip out and over the nozzle nose and/or breaks up at the nozzle 
nose. All the time that these two respective structures are moving downstream, one can observe new structures 
arriving from upstream to replace them. Clearly, the process has a lot of noise and/or velocity fluctuations, but some 
general trends are evident. 
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Figure 15: RSRM Center and Aft Inhibitor Region Vorticity at 100-Second Burn Time 
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Figure 16 shows six frames for the internal motor pressure that correspond to the same times as the vorticity Fig. 
18. These frames were constructed by subtracting out the mean pressure drop as a function of axial distance down 
the length of the motor. So the pressure shown is the fluctuation about the mean pressure for a given axial distance. 
The mean pressure drop down the RSRM bore at 100 seconds is about 7 psi, which is larger than the fluctuations. 
Close inspection of Fig. 16 reveals small low-pressure eddy centers (small green dots) traveling downstream of the 
center inhibitor. Likewise, a low-pressure area (large blue area) appears in the aft segment, but never in the head-end 
in all the figures. This is because the opportunity for large vorticity roll-up exists in the aft-end, but not in the head- 
end. 



Frame No. 1 Run Time = 2.4835 seconds 



Frame No. 2 Run Time = 2.4955 seconds 
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Figure 16: RSRM Center and Aft Inhibitor Region Pressure Contour at 100-Second Burn Time 
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IV. Conclusion 


The complex interactions between internal motor pressure oscillations resulting from vortex shedding, the 
motor’s internal acoustic modes, and the motor’s structural vibration modes were assessed for RSRM and the ETM- 
3 motors. Two approaches were applied to this problem as previously described. Figure 17 presents a comparison of 
the pressure amplitude results from both approaches as compared to measured RSRM data at the 100-second bum 
time. PNCAC001 measured the full pressure (quasi-steady-state and dynamic signal) and PNCAC024 measured 
only the dynamic signal. Results compared reasonably well considering the numerous assumptions involved in the 
two approaches. Figure 18 compares results for the ETM-3 motor at 100 seconds. 
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Figure 17: Comparison of RSRM Head-End Pressure Results at 100-Second Burn Time 
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Figure 18: Comparison of ETM-3 Head-End Pressure Results at 100-Second Burn Time 
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